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Composite pulse sequences designed for nuclear magnetic resonance experiments are currently 
being applied in many quantum information processing technologies. We present an analysis of a 
family of composite pulse sequences used to address systematic pulse-length errors in the execution 
of quantum gates. It has been demonstrated by Cummins et al. in @] that for this family of 
composite pulse sequences, the fidelity of the resulting unitary operation compared with the ideal 
unitary operation is 1 — Ce 6 , where e is the fractional error in the length of the pulse. We derive 
an exact expression for the 6 th order coefficient, C, and from this deduce conditions under which 
this 6 th order dependance is observed. We also present new pulse sequences which achieve the same 
fidelity. 

PACS numbers: 03.67.-a 



Systematic errors are always present in any experimen- 
tal set-up. They are best dealt with by stripping down 
the experiment and finding them one by one. Even after 
doing this, there is still the possibility of not eliminat- 
ing them all. Quantum computers operate universally 
with single-qubit rotations and a controlled NOT gate 
between two qubits [jj. Recently, Cummins et al. Q 
have addressed the issue of systematic errors in single- 
qubit gates due to pulse-length and off-resonance effects. 
Off-resonance errors result in a rotation around an axis 
tilted with respect to the desired rotation axis. Pulse- 
length errors result in a rotation through an angle which 
falls short of, or goes beyond, the desired angle of rota- 
tion, due to an error in the timing of the pulse. Com- 
posite pulse techniques have received much attention as 
a means to correct systematic errors 0, IE IE 13 with 
a significant recent development. First developed in 
the NMR setting, these techniques can be very useful in 
reducing the effects of systematic errors in a wide range 
of quantum computer proposals. They have been identi- 
fied as such in numerous quantum informationprocessing 
technologies such as trapped-ion technologies [USE!) m 
rare-earth doped crystal technology |lllll2| . in supercon- 



ducting technologies 0, Q and in solid state quantum 
information processing technologies [TEl IT^ |. The BB1 
composite pulse sequence, first developed by Wimperis 
0, deals with pulse-length errors in a remarkably effi- 
cient manner. It is shown in that the composite pulse 
sequence has a fidelity of I — Ce 6 when compared to the 
exact qubit rotation. In this article we examine the origin 
of this 6 th order dependence on the error in the fidelity 
of the gate, and deduce some constraints on the possible 
composite pulse sequences one may use. We also suggest 
some other pulse sequences whose fidelities display this 
6 th order dependence. 
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A general single-qubit rotation around an axis in the 
XY-plane of the Bloch sphere can be written as, 



R(6,a) 



exp 



cos a 



Y sin a) 



(1) 



where 9 is the angle through which the qubit is rotated 
and a is the angle of the axis in the XY plane of the 
Bloch sphere, a = being the X axis. We now consider 
the effects of a systematic error, such as a pulse-length 
type error. Using a superscript to denote pulse sequences 
suffering this type of error, where all rotation angles are 
altered by the fractional error e, the above one single- 
qubit rotation becomes 

R e {9,a) = exp ( -i -^ + ^ (X cos a + Y sin a) ) ,(2) 



where e is the fractional error. In order to compare the 
error-prone unitary transformation (V = R e {9,a)) with 
the error-free unitary transformation (U = R{9,a)) 1 the 
following fidelity definition is used, 
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The composite pulse sequence as presented by Wimperis 
takes the form, 

W e {4>u4>2) = i? e K0i)i? e (2^, ( /» 2 )i? e (7r,0 1 ), (4) 

so that when e = 0, W is simply the identity. When this 
pulse sequence is carried out before or after the desired 
single-qubit rotation, the fidelity of the resulting compos- 
ite pulse sequence is I — Ce 6 . In fact, W c (0i,</>2) may 
be placed at any point during a rotation around a given 
axis, i.e., 

R e {a9, a) <-> W £ ((f>i, <fo) <-» i? e ((l - a)9, a), 
with a G [0, 1], and the same fidelity is observed. 
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In order to show where this comes from, we need to ex- 
amine the definition of the fidelity. First we let, 

A = X cos fa + Y sin fa 
B = X cos fa + Y sin fa 
C = X cos a + Y sin a 

be three axes in the XY-plane of the Bloch sphere. C is 
the axis around which we wish to perform the rotation 
and A and B are the two axes in W e (fa, fa). The fidelity 
is, 



Letting P(e) = sbc/i(|; Q, jsbch(e; R, S)), we find that, 

P(e) = e(±Q + R + S) + e 3 (~[R + 2S,[R,S\] 

--L[Q,[Q;2(R+S)]-^[R + S,[Q,R + S]} 
+0(e 5 ). 

Given that all the operators in P{e) are proportional to 
the Pauli operators, Tr(P(e)) — and therefore, 



T = l\Tr{e- i! ^ c WZ{fa,fa)e iiC )\ 

= \\Tr(e-^ c e-^ c W e (fa,fa)e^ c )\ 

= \\Tr(e-^ c W*(fa,fa)e-^ c )\. 

The last step is due to the trace property of invariancc 
under cyclic permutations. W e (fa, fa) can be simplified 
to, 



Tr(e p W) = Tr(l)+Tr(-(P(e)) 2 ) 



= 1 



:Tr((P(e)) 2 ) 



For the next step we let the first derivative of the total 
pulse sequence with respect to the error, e, equal the zero 
matrix at e — to find a relation between the rotation 
axis, C, and the axes A and B, 

The total pulse sequence is given by, 



W e 



, , , -i *V+ c ) A (— i(l+eW.B) 

n, fa) = e 2 e y y ' >e 2 

= e - i *r A {-iA)e- ienE '{-I)(~iA) e - 1 ^ 
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since A 2 = I = B 2 . ABA = X cos(2 fa- fa) +Y $,in{2 fa- 
fa) is another axis in the XY --plane of the Bloch sphere. 
The fidelity now takes on a more symmetric look when 
written as, 



T = L\Tr(e~ l ^ c e~ l ^ A e- l£nABA e-' l ^ A e-^ c )\ 



This form for the fidelity turns out to be the reason the 
BB1 sequence performs so well. The symmetric Baker- 
Campbell-Hausdorff formula is stated and derived in Q| 

as, 



e \tR^,S e \tR _ e sbch(t-R,S) 



(5) 



with, 



sbch{t;R,S) = t(R + S) - —t 3 [R + 2S,[R, S]} + 0(t 5 ) 

If we define Q = -i9C, R = -iirA and S = -iirABA 
to simplify the notation, we can reduce the expression in 
the trace above by twice applying the symmetric BCH 
formula, yielding, 

T = ^\Tr(e-^ c e~ l ^ A e~ tenABA e- t ^ A e-^ c )\ 



\Tr(e^e^ R e tS ei R e^)\ 



= t\Tr(ei Q e^ sbch{t '' R ' S)) e^ Q )\ 



■|Tr(e 



sf)c/i(f ;Qlsbch{e;R,S))\ 
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so that, 



e iipl Qe iipl Re(1+t) S e l±plR 



dBB(e) 



dc 



{R + S) = 0. 



-=o 



This allows us to simplify the expression for P(c) 
P(e) 



e\-h,R + 2S, [R,S]] + ^{Q,{Q,Q]} 



1 

96 
1 

24 



[Q,[Q,Q]]) + 0(e 5 ) 

[R + 2S, [R, S]]e 3 + 0(e 5 ). 



Next we expand the commutator [R + 25, [R, S]] to get, 

{\R + 2S,[R,S}}) 2 = -tt 6 ([A + 2ABA, [A, ABA}}) 2 
= -tt 6 (40I + 8(AB + BA) 
-20(ABAB + BAB A) 
-8(ABABAB + BAB ABA)) 
= -7T 6 [40 + 16cos(0 2 - fa) 

- 40 cos 2(02 - fa) 

- 16cos3(0 2 - faW 

since AB = cxp(i(fa — fa)Z). 
Therefore, we can write the fidelity as, 



T = 1 
= 1 



1 



2304 



144 



e 6 Tr({[R + 2S, [R, S}}) 2 ) + 0(e 8 ) 

[5 + 2cos(0 2 - fa) - 5 cos 2(02 - fa) 
-2cos3(fa-fa)] + 0(e s ) 
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Finally, by writing A, B and C in terms of the Pauli 
operators, X and Y, the condition + R + S = 
implies, 



cos 0i + cos(20i — 2 ) = ——cos a (6) 

Q 

sin 0i + sin(20i — 02) = sin a, (7) 

In 

and these two equations lead to the following rule for 
choosing 0i and 02 given 9 and a, 



h = a — - arccos I — - - 1 
2 \ 87r^ 

!>2 = 30i — 2a 



(8) 
(9) 



To give a concrete example, we choose a = 0, 9 = n 
corresponding to a 180° rotation about the X axis of the 
Bloch sphere. This gives us, 



1 



30i 



and a fidelity, 



T = 1- 



1024 



which agrees with the findings in 

It can be readily seen that if n copies of W e (0i , 02 ) are 
carried out one after the other, then the condition that 

dBB I 
de I e=0 



0. 



1 



-Q + n(R+S) = 



and also P(e) = e(±Q + n(R + S)) + 0(e 3 ). This defines 
the general Wn pulses introduced in [2j as, 



a arccos 

2 

30i - 2a 



8n 2 7r 2 



and shows that they all perform with a fidelity 1 — 0(e e ). 
We note that it is only necessary to show the sequence 
works for one angle, say a = 0, with other axes accounted 
for by phase-shifting each pulse by some angle a. 



Other 6 th Order Pulse Sequences 

We now look at whether there are other types of pulse 
sequences similar to those above which can achieve the 
same fidelity. The most general 3-pulse sequence is, 



W e 



,0 2 ,0 3 ) = i? £ (C,0i)i? e (77,0 2 )i? £ (7,03). 



As before we define axes A, B and C in the XY-plane of 
the Bloch sphere by the angles 0i , 02 and a respectively. 
The first condition to satisfy is that the pulse sequence 
must be symmetric so that it can be reduced using the 
symmetric BCH formula, thus keeping the fidelity for the 
total pulse sequence 6 th order in e. This means we need 
C = 7 and 0i = 03. Also to satisfy the condition that 
Vy e (0i, 02) = I, we require, 

2-f + ri = 4mir,m = 1,2, ... (10) 

We are then left with sequences of the form, 



W e 



,0 2 ) = i? e ( 7 ,0i)i? £ (2(2m^- 7 ),02)i? £ ( 7 ,0i) 



Finally, we let the first derivative of the total pulse se- 
quence with respect to e equal zero at e = and obtain 
the constraint, 



= -C + jA + (2m7r - 1 )e~ i 'i A Be 1 ^' 



sin 7 sm(0 2 



= o, 



giving 7 = p7r, rj — 2q~K with p + q = 2m and 
p,q = 1,2,... The most general 3-pulse sequences are 
R e (pir, 0i ) R e (2qir, 2 ) R e (pir, 4>i) with 0i, 2 determined 
from, 



2tt 



P A + qA p BA p = 0. 



For even p, 



pcosi 



q sin <p 2 



q cos <p 2 



2tt 



2tt 



sin(a) 
cos(a) 



meaning that either 0i 



arcsin(| sin(0 2 



a)) or 



02 = a + arcsin(^ sin(a — 0i)). Hence p/q < 1 since 
~ 1 < g sin(a — 0i) < 1 Va. Similarily, q/p < 1, i.e. 
q = p and 0i + 02 = 2a. The same is true for odd p 
except 30i — 02 = 2a. 

We are left with the general 3-pulse sequences, 

W4(0i,0 2 ) = R e (mn,^ 1 )R'(2m 1 TT,(j )2 )R e {mTT,(j )1 ) 

The Wn family is obtained from repeating the 
Fyi(0i,0 2 ) sequence n times. When m = 2, the "pass- 



band" Wimpcris sequence [3J, |6j is recovered. Moreover, 
these are the only 3-pulse sequences which achieve this 
Qth or( j er fidelity. The first three of these composite se- 
quences are plotted in Fig. ^ In other 3-pulse sequences, 
the first order term in the fidelity always disappears due 
to the fact that by collapsing the entire pulse sequence 
using the BCH formula to, 



BB(e) 



«g(l + e) 

e — W e 



2) 
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the fidelity is T = 1 + \Tr((P\ (e)) 2 ) + • • • and the leading 
term in P\{e) is 0(e). 

So there are no pulse sequences constructed from 3 
pulses which achieve this 6 th order dependence for the 
fidelity for the resulting rotation other than the above 
family of pulse sequences. There is, however, the option 
of creating a 5-pulse sequence by introducing a third axis. 
The general form for such a sequence is, 

i? e (C, <k) R f (v, h) i? e (7, h) R e (», R e (v, 5 ). 

We define the axes A,B,C and D in the AF-plane of the 
Bloch sphere by the angles 0i , 2 , 03 an d Oi respectively. 
D is now the axis around which we wish to rotate by 9. 
Again, we will require £ = v,r) = \i and 0i = 05, 02 = 04 
in order to keep symmetry in the sequence and hence the 
6 th order fidelity dependence. In order that the pulse 
sequence is the identity when e = 0, we need to satisfy, 
2Q + 2?7 + 7 = 4m7r, m — 1, 2, ... The first derivative of 
the total pulse sequence is zero at e = when £ = pn, 
rj = qn for positive integers p and q. We arrive at general 
5-pulse sequences of the form, 

W e = 

i? e (p7r, t ) i? £ (g7r, fa) R e (2rir, 3 ) R £ (qir, 2 ) R £ (pn, 0i). 

where p + q + r = 2m and 0i , 02 and 03 are determined 
from, 

—D + pA + qA p BA p + rA p B q CB q A p = (11) 
2tt 

One solution to find 5-pulse sequences is at p = q = r, 
analogously to the 3-pulse sequence case. In this case 
p = and, as it must remain an integer, m must be a 
multiple of 3 with now p = 2, 4, 6.... Eqn. (11) is satisfied 
for p = 2 when 0i = 0, 02 = arccos( e ~^ 7r ) and 03 = —02 
for a rotation around the -X-axis. However, while the 



fidelity displays a 6 order dependence on the fractional 
error , the coefficient of the leading term, shown in Table 
HI is so much larger that the sequence is only better than 
the error-prone pulse for small values of e (~ 0.2). The 
situation does not improve for higher values of p and so 
these sequences are of no real practical use. 

Another 5-pulse sequence which achieves the same 6 th 
order dependence for the fidelity is found by setting 
p = 1, q = 2, r = 1. Eqn. (11) is now satisfied when 
0i = arccos( 6> ^ 7r ), 02 = 20i, 03 = 30i for a rotation 
around the -A-axis (a = ir). As before, other axes may 
be accounted for by phase-shifting each pulse by the ap- 
propriate angle. The fidelity of this sequence is much 
better than the previous 5-pulse sequence and is quite 
close to that of the PB1 sequence as seen in Fig. ^ 
Other sequences can be constructed by varying p, q and 
r. 



Conclusion 

We have presented an analysis of the composite pulse 
sequences presented by Jones et al. 0, to combat sys- 
tematic pulse-length errors in single-qubit rotations. We 
have derived an explicit form for the fidelity and shown 
how it is possible to set up other 3-pulse sequences which 
achieve the same order error dependence for the fidelity. 
We have shown that there are also 5-pulse sequences 
which do achieve the 6 th order dependency of the fidelity 
on the error. 
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59.1 




190.6 


wi 


283.4 


wi 22 


877.8 



TABLE I: The coefficients C in the fidelity expansion F — 
1 — Ce 6 for six composite pulse sequences which compensate 
for an error-prone 7r-pulse around the -X-axis. 




FIG. 1: Fidelity of composite pulse sequences for (a) the 
3-pulse sequences (solid) and (b) the 5-pulse sequences 
(dashed), with coefficients for the 6 th order term of the fi- 
delity given in Table 1. The fidelity of the single error- prone 
pulse is also shown (dotted). 



